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Abstract 

The  NLO  properties  in  terms  of  static  and  frcquency'dcpcndcnt  hyperpolarizabilities  of  Si  (OH ),  and  Ge(OH)4  clusters  are 
determined  with  an  aim  to  understand  the  microscopic  NLO  behavior  of  silica.  The  calculated  results  reveal  that  the  cluster 
configuration  induced  by  the  presence  of  a  defect  plays  a  significant  role  in  determining  the  non-linear  response  of  these  clusters. 
In  addition,  instant  ionization  of  the  cluster  introduces  a  noticeable  increase  in  the  value  of  second  byperpolarizability,  y. 


Nonlinear  optics  is  currently  at  the  forefront  of  re¬ 
search  because  it  is  the  key  technology  for  optical 
communication,  signal  processing  and  computing 
[  1  ].  Photonic  switches  made  out  of  nonlinear  optica] 
(NLO)  materials  have  potential  advantages  of  wide¬ 
band  frequency  domain,  immunity  to  electromag¬ 
netic  induction,  and  high-speed  operation  over  elec¬ 
trical  switcltes.  Glassy  materials  are  one  of  the  prime 
candidates  for  photonic  switches  due  to  their  case  of 
fabrication,  high  transparency  and  high  chemical  and 
thermal  stability  [2], 

Nonlinear  optical  properties  of  glasses  are  known 
to  be  correlated  with  the  third-order  susceptibility  due 
to  the  absence  of  ordered  dipoles.  However,  Kester 
ct  al.  [3]  have  observed  frequency  doubling  in  Ge- 
doped  silica  suggesting  that  even  glasses  have  fre¬ 
quency  doubling  capability.  Similar  observations  have 
been  made  in  bulk  and  thin  films  of  silica  grown  on 
semiconductor  substrates  [4]. 

The  present  work  attempts  to  provide  a  basis  to 
understand  the  origin  of  such  NLO  eficcts  in  pure  and 
impurity-doped  silica  at  the  molecular  level.  Our  ef¬ 
forts  will  be  directed  to  hyperpolarizability  calcula¬ 
tions  of  the  molecular  units,  namely  Si(OH)«  and 


Ge(OH)4  clusters  (Fig.  1)  in  the  Hartrce-Fock 
framework.  Each  oxygen  atom  at  the  cluster  bound¬ 
ary  is  attached  by  a  hydrogen  atom  in  such  a  way  that 
they  simulate  0-Si  bond  directions.  The  optimum 
configurations  of  these  nine-atom  clusters  are  ob¬ 
tained  using  a  general  atomic  and  molecular  elec¬ 
tronic  siructiirc  program  [5],  GAMESS. 

Tlie  NLO  properties  have  been  shown  lo  be  very 
sensilive  lo  the  choice  of  basis  set.  For  the  case  of 


Fig.  I.  Geometrical  configuration  of  the  X(OH)4  cluster  where 
X  is  Si  or  Gc. 
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ethylene,  a  sufTiciently  flexible  valence  Gaussian  ba¬ 
sis  set,  such  as  a  3-21 G  or  6-3 IG,  augmented  by  dif¬ 
fuse  p  and  d  functions  on  C  atoms  is  found  to  be  nec¬ 
essary  for  a  satisfactory  description  of  various  tensor 
components  of  the  static  polarizability  and  second 
hypcrpolarizability  [6].  ForSi(OH)4and Gc(On)4 
clusters,  we  use  double-zeta  valence  (dzv)  basis  sets, 
with  a  diffuse  d  polarization  function  to  describe 
cluster  atoms.  The  basis  set  consists  of  6s  and  4p  and 
Id  Gaussian  functions  for  Si  [7],  6s,  4p  and  2d  for 
Ge  [8], 3s,  2p  and  IdforC  [7]  and  2sforH  [7], 

For  the  dipole  moment  of  a  molecule  interacting 
with  a  static  electric  field  (F),  we  can  write 

III  —  Itif)  "I"  ^'/  ^i 

+  i)'i(*/FyF* .  (I) 

where  /(,o  is  the  permanent  dipole  moment  and 
flijh  and  yijki  are  tensor  elements  of  the  linear  polariz¬ 
ability,  and  first  and  second  hyperpolarizabilitics.  re¬ 
spectively,  of  the  molecule  [9].  Likewise,  one  can 
write  an  expression  for  dynamic  fields  relating/?  and 
)'  to  experiments  such  as  second  harmonic  genera¬ 
tion,  optical  rectification,  and  Kerr  effect  [  10].  In  Eq. 
( I ),  the  subscripts  i,  j,  k,  and  /  represent  Cartesian 
axes  and  the  summation  of  repeated  indices  is 
assumed. 

First,  we  obtain  the  potential  energy  surfaces  of 
both  neutral  (geom  A,  q=0)  and  ionized  states 
(geom  B,  +  1 )  ofthe  X(OH)4cluster  (X  being  Si 


Fig.  2.  Sketch  of  the  potential  energy  surface  of  the  neutral  and 
ionized  X(OH  >4  cluster  where  X  is  Si  or  Ge. 

orGc)  and  then  calculate  their  hyperpolarizabilities 
at  the  optimum  configuration.  The  ^=-t-l  charge 
state  is  expected  to  mimic  the  defect-induced  effect 
on  the  NLO  properties.  We  also  use  the  vertical  ap¬ 
proximation  for  calculations  of  the  single  point  A' 
(geom  A,  ^  =  +  1 )  and  B'  (geom  B,  ^=0)  as  shown 
in  Fig.  2. 

Table  1  shows  the  geometric  parameters,  namely 
the  bond  length  and  the  bond  angles  of  the  optimized 
configuration  obtained  with  dzv  and  dzvP  (includ¬ 
ing  polarization  function,  d)  basis  sets.  The  results 
obtained  using  both  the  basis  sets  show  similar  trends. 
For  example,  the  substitution  of  Ge  for  Si  increases 
the  bond  length.  Mulliken  population  analysis  pre¬ 
dicts  Ge  to  be  less  electronegative  than  Si.  For  the 
ionized  cluster  (q=  +  1 ),  the  optimized  configura¬ 
tion  shows  the  hole  localizing  on  one  of  the  neighbor¬ 
ing  oxygens  resulting  in  an  increase  of  its  bond  length. 


Table  I 

Total  energy  and  configurational  parameters  of  neutral  and  ionized  Si(OII  )4  and  Gc(OI  1  >4  clusters 


/t(X-0(i),i  =  2,  5) 

(A) 

/l(0(2)-X-0(i).  1  =  3,  5) 
(deg) 

Total  energy 
(hartree) 

Si(OH)4 

<7=0 

dzv 

1.652,  1.632.  1.652,  1.652 

107.3.  113.8.  107.3 

-590.8126 

dzvP 

1.628,  1.628,  1.628,  1.628 

106.5.  115.7,  106.5 

-590.9846 

6-31G  (llj 

1.665,  1.665,  1.673.  1.673 

115.5.  107.8.  106.2 

-590.7039 

<7=  +  l 

dzv 

1.594,  1.600.  1.584.  1.943 

117.8,  119.4.95.6 

-590.4592 

dzvP 

1.578,  1.582,  1.575,  1.900 

118.2.  119.8.95.0 

-590.5971 

6-3IG  [It] 

1.602,  1.918,  1.615,  1.608 

102.8.  116.9.  120.2 

-590.3531 

Ge(OH)4 

II 

0 

dzv 

1.726,  1.726,  1.726,  1.726 

107.6,  113.3,  107.6 

-2375.0059 

dzvP 

1.734,  1.734,  1.734.  1.734 

106.3.  115.9,  106.3 

-2375.3107 

9=  +  l 

dzv 

1.655,  1.653,  1.654,  1.953 

117.7.  118.6,  96.4 

-2374.6614 

dzvP 

1.694,  1.691,  1.692,  1.999 

II6.I.  II6.I.9I.6 

-2374.9574 
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For  comparison,  we  include  the  results  of  Rao  et  al. 
( 1 1  ]  using  the  6-3 1 G  basis  set  for  the  Si  ( OH  )<  clus¬ 
ters  in  Table  1. 

In  order  to  judge  the  character  of  predicted  config¬ 
uration  on  the  potential  energy  hypersurface,  we  cal¬ 
culated  the  frequency  of  the  neutral  Si(OH)4  cluster 
(i.e.  geomA)  in  the  harmonic  approximation.  The 
calculated  frequencies  in  the  rocking,  bending  and 
stretching  modes  come  out  to  be  465,  828  and  1041 
cm"',  showing  a  good  agreement  with  the  corre¬ 
sponding  experimental  values  [12]  of  430,  812  and 
1058  cm"',  respectively. 

Table  2  lists  dipole  moment  and  the  hyperpolariz¬ 
abilities  for  the  Si(OH)4  and  Ge(OH)4  clusters.  As 
we  ionize  the  cluster  (freezing  the  configuration  at 
9=0),  the  P  and  y  show  a  large,  noticeable  increase. 
The  optimized  configuration  (for  <7=  -I- 1 )  retains  the 
somewhat  larger  P  value.  However,  when  we  make 
the  cluster  charge  neutral,  both  the  dipole  moment 
and  ^  show  an  increase  in  their  values.  Similar  obser¬ 
vations  have  been  made  for  the  Ge(OH)4  cluster. 
Note  that  for  q=  + 1  (at  ^=0  configuration)  a  large 
increase  of  y  is  predicted. 

The  results  therefore  reveal  that  the  cluster  config¬ 
uration  plays  a  key  role  in  predicting  a  large  p  value 
related  to  second-harmonic  generation  or  Kerr  effect. 
A  dependence  of P  on  the  cluster  dipole  moment  re¬ 
sulting  from  the  distorted  configuration  is  also  shown 
here.  On  the  other  hand,  instant  ionization  of  the 
cluster  leads  to  a  significant  increase  in  y  that  is  re¬ 
lated  to  DC-SHG  and  THG  experiments. 

Tabte  2 

Dipole  moment  and  polarizabilities  (in  atomic  units)  of  neutral 
and  ionized  Si(OH)4  and  Ge(OH)4  clusters  at  diflerent  geome¬ 
tries  shown  in  Fig.  2 


Dipole  moment 

y 

Si(OH)4 

geom  A  (9=0) 

0.00 

0.00 

945.8 

A'(9=-H) 

0.00 

5.84 

22420.8 

B  (9= -1-1,  opt) 

0.24 

31.81 

3008.9 

B'(9=0) 

1.67 

63.6 

1455.1  ■ 

Ge(OH)4 

geomA  (9=0) 

0.00 

0.23 

1248.1 

A'(9=  +  1) 

0.00 

34.5 

660665.5 

B  (9=  4- 1,  opt) 

0.59 

1.42 

1690.9 

B'  (9=0) 

1.36 

9O.C1 

2226.3 

Table  3 

1'inie-dcpcndcnt  tlarlrcc-Fock  results  for  the  Si  (OH  >4  cluster  al 
IU64  nm  in  atomic  units 


Geom  A 

Geom  B' 

SHG 

0.1 1 

63.8 

Kerr  effect 

0.1 1 

61.5 

DC-SHG 

1185.7 

1713.4 

THG 

1254.5 

1848.7 

It  has  been  suggested  that  the  dispersion  of  NLO 
effects  is  of  great  importance  for  the  application  of 
conjugated  polymers  to  photoionic  studies.  But  this 
is  not  found  to  be  the  case  for  both  the  Si  (OH  >4  and 
Ge(OH  >4  clusters.  We  give  the  results  of  time-depen¬ 
dent  Hartree-Fock  (TDHF)  calculations  using  the 
program  package  ACES  II  [13]  al  1064  nm  in  Table 
3.  Although  the  frequency-dependent  values  arc  pre¬ 
dicted  to  be  about  the  same  as  the  static  values,  they 
do  confirm  the  role  of  configuration  for  a  larger  P 
value.  Our  results  therefore  suggest  that  the  cluster 
configuration  induced  by  the  presence  of  a  defect 
plays  a  role  for  the  NLO  effect  at  the  molecular  level 
in  the  Si(OH)4and  Gc(OH)4  clusters. 

Helpful  discussions  with  Dr.  John  Kester  are  ac¬ 
knowledged.  We  would  like  to  thank  the  Air  Force 
Materials  Command,  the  Air  Force  Office  of  Scien¬ 
tific  Research,  and  Division  of  Chemical  Sciences  of 
the  Seiler  Research  Laboratories  at  the  USAF  Acad¬ 
emy  for  sponsorship  of  this  research. 

References 

1 1 )  J.D.  Mackenzie,  C.Y.  Li  and  Y.  Xu,  Chem.  Express  6  ( 1 99 1 ) 
903. 

[2)  E.J.  Friebele  and  D.L.  Griscom,  Mater.  Res.  Symp.  Proc. 
61  (1986) 319. 

1 3  J  J.  Kester  el  al.,  private  communication. 

(4)R.A.  Myers,  N.  Mukherjee  and  S.R.  Brueck,  Nonlinear 
Optics:  Topical  Meeting,  August  17-21,  1992,  Hawaii, 
abstract  MA3. 

|S]  M.W.  Schmidt,  K.K.  Baldridge,  J.A.  Boatz,  J.H.  Jensen,  S. 
Koseki,  M.S.  Gordon,  K.A.  Nguyen,  T.L.  Windus  and  S.T. 
Elbcrl,  QCPE  Bull.  10(1990)  52. 

16]  S.P.  Kama,  G.B.  Talapatra,  W.M.K.P.  Wijekoon  and  P.N. 
Prasad,  Phys.  Rev.  A  45  ( 1992)  2763. 

[7]  T.H.  Dunning  Jr.  and  P.J.  Hay,  in:  Methods  of  electronic 
structure  theory,  ed.  H.F.  Schaefer  III  (Plenum  Press,  New 
York,  1977)  p.  1. 


340 


R.  Pandey  et  al.  /  Chemical  Physics  Letters  220  (1994)  S37-340 


{8]R.C.  Binning  Jr.  and  L.A.  Curtiss,  J.  Compul.  Cheni.  II 
(1990)  1206. 

(9]  C.  Flytzanis,  in:  Theory  of  nonlinear  optical  susceptibilities, 
quantum  electronics,  eds.  H.  Rabin  and  C.L.  1'ang 
(Academic  Press,  New  York,  1 975). 

( lO)  D.S.  Chemla  and  J.  Zyss,  eds..  Nonlinear  optical  properties 
of  organic  molecules  and  crystals,  Vols.  I  and  II  (Academic 
Press,  New  York,  1987). 


(Ill  P.S.  Rao,  R.J.  McEachern  and  J.A.  Weil,  J.  Comput.  Chem. 
I2 (I99I)  254. 

( 1 2 )  H.M.  Dauplaise,  K.  Vaccaro,  B.R.  Bennett  and  J.P.  Lorenzo, 
J.  Elecirochem.  Soc.  1 39  (1 992)  1 684. 

|I3)  H.Sekino  and  R.J.  Bartlell,  J.  Chem.  Phys.  85  (1986)  976; 
H.  Sekinoand  R.J.  Bartlett.  J.  Chem.  Phys.  94  (1 99 1 )  3665; 
J.F.  Stanton,  J.  Gauss,  J.D.  Walts,  W.J.  Lauderdale  and  R.J. 
Bartlett,  University  of  Florida  (199I). 


